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Abstract 

Intransitivity is a property of connected, oriented graphs representing species interactions that may drive their coexistence 
even in the presence of competition, the standard example being the three species Rock-Paper-Scissors game. We consider 
here a generalization with four species, the minimum number of species allowing other interactions beyond the single 
loop (one predator, one prey). We show that, contrary to the mean field prediction, on a square lattice the model 
presents a transition, as the parameter setting the rate at which one species invades another changes, from a coexistence 
to a state in which one species gets extinct. Such a dependence on the invasion rates shows that the interaction graph 
structure alone is not enough to predict the outcome of such models. In addition, different invasion rates permit to tune 
the level of transitiveness, indicating that for the coexistence of all species to persist, there must be a minimum amount 
of intransitivity. 
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1. Introduction 

Cyclic competition [ll-Q among a population of S 
species (or S different traits within a species) may occur 
when the trophic network presents loops, for which several 
examples exist: mating lizards competing bacteria [3- 
Q , coral reef environments Q , competing grasses Top 12| , 
etc. The simplest and most studied case corresponds to 
the Rock-Scissors-Paper (RSP) game, with 5 = 3, in 
which each strategy dominates the next one, in a cyclic 
way [3, 14 1 . These interactions, or food chain, are thus 
given by a three vertices, single looped oriented graph. 
Since there is no perfect ranking of the species, the system 
is fully intransitive. A direct generalization [15l4l9j is to 
consider S > 3 competitors whose interactions also follow 
an oriented ring, — >■ 1 — > ... — > S — 1 — » 0. For the spe- 
cific case of S = 4 17H24|. the minimum value for which 



neutral pairs may exist, those non interacting alliances 
help prevent invasions. Such defensive alliances may also 
appear between non mutually neutral species (cyclic al- 
liances) when the interaction graph has more than a single 
loop 20- 22, 24 - 37 1 . Random and non regular food webs 
have also been considered [38|-|40|. 

Such models, with simplified competing interactions and 
food webs, do not claim quantitative predictions, but at- 
tempt instead to unveil the universal behavior that results 
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from the direct competition between species. The inter- 
actions are coarse grained in the sense that the ultimate 
mechanism (dispute for space, resources, mating partners, 
etc) and its non all- nothing nature (e.g., dependence on 
size, age, distance and other contingent factors) are aver- 
aged out and replaced by a simple, probabilistic interac- 
tion. Such interactions may depend on space, time, be a 
characteristic of the two species involved, etc., what intro- 
duces heterogeneities in the system 17 - 1|| 2g, 41 - 4fj| . In 
turn, this gives rise to hierarchical alliances and diverse 
levels of intransitivity. Anomalous, negative responses 
may occur in this case, an example being the "survival 
of the weakest" principle, observed for 5 = 3 4l|, |47[ (and 
its generalization for S > 3 [HI, [l3|), in which a species 
density may increase after its invasion capability has been 
decreased. Real systems, with their more complex trophic 
networks, may even have more complex responses to vari- 
ations in the invasion rates and, consequently, predicting 
their behavior in such a situation will be far more difficult. 

Spatial correlations may exist when the range of inter- 
action is limited but play no role when the interactions 
are spatially unconstrained (fully mixed case), and simple 
mean field approximations are expected to produce reason- 
able results for sufficiently large systems in such a case. 
Nonetheless, stochastic fluctuations are expected to be- 
come important for finite size populations and even drive 
the system towards one of its absorbing states, in which 
one or more species become extinct, decreasing the diver- 
sity. 

Intransitivity is considered a key mechanism for diver- 
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sity sustaining in the presence of competition. Thus, im- 
portant questions arise on the effects of tuning the transi- 
tivity by changing the invasion probabilities. For example, 
does the diversity suddenly decrease once the system is no 
longer fully intransitive? Can diversity be predicted solely 
based on the structure of the interaction graph? How does 
the system respond to changes on the interaction param- 
eters of a complex trophic network? 

To answer these questions, we start with a fully intransi- 
tive ring of four species (S = 4) competing with the same 
unity invasion rate. All four species have similar roles, 
with one prey and one predator each. This symmetry is 
broken when the interaction graph is turned into a fully 
connected graph, with two diagonal interactions having a 
rate x of invasion, as shown in Fig. [1] This introduces 
some hierarchy in the system: the top species, and 1, 
have two preys each (and one predator), while the bottom 
ones, 2 and 3, have two predators (and one prey). The 
arrows indicate the direction in which the invasion occurs 
and the corresponding rate: around the original ring, in- 
vasions occur with unitary rate while along the diagonals, 
this probability is \- Species 2 and 3 have only one prey 
each, but once they encounter their prey, they always sub- 
jugate them. On the other hand, species and 1 have 
two preys each, but with a smaller than unity success rate 
and are the weakest species in this case (see Section 2] for 
the detailed discussion). When \ = 1 we recover the case 
considered in Ref. [35[ (see also Ref. (2ll|V 





Figure 1: Species interactions represented by an oriented graph 
(right), the arrows indicating the invasion direction and 1 and 
< x ^ 1 are the corresponding rates. The possible configurations, 
(pOi Pi i P2i Ps), can be represented as points in a 3-simplex (left). 
Upon extinction of a single species the orbits are constrained to one 
of the faces of the tetrahedron, a 2-simplex. After a second extinc- 
tion, the configurations are constrained to the line joining the two 
surviving species (1-simplex). If eventually a single species remains, 
this state is represented by the corresponding vertex. 



The paper is organized as follows. Next section dis- 
cusses the mean field approach for the fully mixed version 
of the model, in particular, the stable fixed points both for 
X = and v^O. Then, we present the results for the spa- 
tially structured system, with emphasis on the long term 
persistence of the coexistence state. Finally, we discuss 
the similarities and discrepancies of both approaches and 
present our conclusions. 



2. Analytical Results 

When spatial correlations are neglected and individuals 
have the same probability to interact with all others, irre- 
spective of their distance, one may attempt a mean field 
description. Let pi be the density of species i (obviously 
Si Pi ~ !)■ Time variations in the densities may only oc- 
cur due to interactions between different species, in which 
the stronger one will invade the weaker with rate 1 or x- 
The mean field equations depend only on the frequency of 
such encounters and read: 
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where each element of the interaction matrix, Jy, is the 
rate with which species i invades j. A negative 7jj means 
that the invasion direction is reversed. The matrix I is 
given by 



(2) 



These equations present several equilibrium points such 
that pi — 0, Mi. The linear stability of these steady states 
is determined by the sign of the real part of the eigenval- 
ues of the Jacobian matrix. If at least one eigenvalue has a 
positive real part, the corresponding fixed point is unsta- 
ble, otherwise it is stable. Furthermore, a stable equilib- 
rium point may be asymptotically attainable when all real 
parts are strictly negative. When there are purely imag- 
inary eigenvalues, the stable equilibrium is neutral and 
never attainable dynamically. 

The fixed points for v = have been discussed by sev- 
eral authors [2, 17 - 19|, [23| . First, there are four absorbing 
states that are heteroclinic points (saddle points) [l|, at 
the vertices of the 3-simplex, in which only one species 
survives: (1,0,0,0), (0,1,0,0), (0,0,1,0) or (0,0,0,1). In 
addition to these, and because species and 2 (or 1 and 
3) are mutually neutral, any point on the line connect- 
ing each pair is also a fixed point, the initial proportion 
between them kept constant: 



(c D , 0,1- Co , 0) 
(0,c o ,0, 1 - c ) , 
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with < c < 1. Lastly, there is a coexistence fixed point 
in the interior of the 3-simplex, for which all densities are 
non zero, 

1 1 
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with < c < 1/2, a particular example being the sym- 
metric state (1/4,1/4,1/4,1/4). This point is stable, but 
not asymptotically stable. In fact, there are two integrals 
of motion: P0P2 and P1/O3. In contrast, for the 5 = 3 game 
there is only one invariant of motion: po/°i/°2- 

With x 7^ 0, the coexistence state Eq. (J5J) is no longer 
a solution of Eq. (jTJ) . Nonetheless, there are two further 
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fixed points in which one species (1 or 2) dies out and 
the remaining three species form a non homogeneous RSP 
game [35j ]: 

Notice that for % = the above fixed points are partic- 
ular cases of Eqs. and (U]), respectively. The first so- 
lution, in which the species 1 becomes extinct, is an un- 
stable fixed point, while the second one, in which species 
2 goes extinct and the remaining three compose a het- 
erogeneous RSP game, is (neutrally) stable. In the limit 
X — > oo, the stable solution Eq. (J7J becomes (1, 0, 0, 0) and 
species dominates. An example, Fig. [21 shows the evo- 
lution from the symmetrical initial state with pi = 1/4, Vi 
and x = 0-5- The system approaches a closed orbit that 
oscillates around (0.2,0.4,0,0.4), Eq. ([7]), after the expo- 
nentially fast extinction of species 2. When P2 = 0, the 
quantity p^p\Pz is an integral of motion [l], |3 [49[. In- 
terestingly, besides the trivial normalization condition, no 
invariant involving all four densities exists for x^0. The 
fixed points are equivalent to the time average of the os- 
cillating densities. Both the period of the oscillations and 
the time that species 2 takes to become extinct diverge 
when x — > since in this limit the existence of the in- 
variants of motion mentioned above precludes the possi- 
bility of an extinction. Indeed, when X = 0, the homo- 
geneous initial condition considered here becomes a fixed 
point with four coexisting species. For the homogeneous 
case, x = 1) one recovers the pi = 1/3, Vi, solution. Notice 
that each species density depends on its prey's invasion 
rate and when we decrease x (species 1 invasion rate over 
3), although one would expect a decrease in its density, 
it is the density of its predator, species 0, instead, that 
decreases. This is known as the "survival of the weakest" 
principle jill l47j| . 

When spatial correlations are important, as when agents 
are placed on a lattice (see Section [3]), the mean field ap- 
proach usually breaks down. One of the simplest ways to 
go beyond the mean field predictions is to use the pair 
approximation (PA) [Hoi ]. Within this approach one con- 
siders the dynamics of pairs of connected sites (instead of 
only one-site quantities as in MF). As the corresponding 
equations depend on triplets of connected sites, the sys- 
tem is closed by choosing an ansatz relating three- and 
two-site quantities. For the PA the ansatz chosen is of 
the form P(123) = P(12)P(23)/P(2), where P(123) is the 
probability of having species 1, 2, and 3 occupying three 
connected sites (2 occupies the central site). P(12) and 
P(23) are similarly defined. 

For the system considered here the PA does not have 
any fixed points with coexistence of the four species. In 
fact, when expressed in term of species densities, the fixed 
points of the PA coincide with those found using MF. One 
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Figure 2: Example of the time evolution for \ = 0.5 within mean 
field, starting with pi(t = 0) = 1/4, Vi. Species 2 rapidly becomes ex- 
tinct and the remaining three species oscillate, out of phase, around 
the respective neutrally stable fixed point (dashed lines) given by 
Eq. (7), (0.2,0.4,0,0.4). The fixed point also corresponds to an av- 
erage of pi over a period, once the stationary state is attained. 



important difference is that the fixed point for which there 
is extinction of species 2 is asymptotically unstable. As 
happens in the case of the RSP game (511 ? ], the model 
has a heteroclinic cycle involving the four species. In addi- 
tion, the diagonal interactions give rise to two new hetero- 
clinic cycles involving species 0, 1 and 3, and 0, 2 and 3. As 
these cycles share some nodes, none of them can be asymp- 
totically stable. This, however, does not mean that they 
cannot dominate the dynamics. Solving the PA equations 
of motion for several different initial conditions shows that 
for long times the system asymptotically approaches the 
3210 cycle. However, for shorter times the system moves 
towards the vicinity of the 310 cycle, and can stay there 
for a long time until it "jumps" to the vicinity of the 3210 
cycle. This can be thought of as a 'competition' between 
the cycles [52| that is eventually won by the cycle 3210 (the 
cycle 023 does not seem to play any role in the dynamics). 
The fact that the density of species 2 falls to extremely 
low levels during the transient implies that in a stochastic 
version of this dynamics the extinction of species 2 would 
happens after a rather short time. The duration of this 
transient is an increasing function of x- Unfortunately, it 
is not possible to obtain the time average of the densi- 
ties of an y o f the species because these quantities do not 
converge [531 ] . The above behavior is illustrated in Fig. [31 
In particular, notice that species 2 resumes to a noticeable 
density in the top panel of Fig. [3] All four densities appear 
on the heteroclinic orbit, one at a time (with the others 
being extremely small) with increasing periods of stasis 
on each (unstable) monoculture [? ? ]. Although all four 
species have non zero densities, the heteroclinic cycle is 
termed "impermanent coexistence" [? ] since they do not 
coexist with finite densities. In the bottom panel, on the 
other hand, the orbit is depicted in the simplex of Fig. [TJ 
Initially (the starting point is the black dot) the orbit ap- 
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Figure 3: (Top panel) Numerical resolution of the pair approximation 
equations for \ = 0.45, starting with pi(t = 0) = 1/4, Vi. Notice 
that f>2 decreases to very small values during the transient, and the 
system approaches the 310 cycle. Later, however, p2 grows once 
again and the system switches to the 3210 cycle. (Bottom panel) 
The same orbit plotted in the simplex. The black dot signals the 
starting point. Although initially the orbit approaches the 310 face 
and the species 2 seems to go extinct, it eventually resumes and the 
full 3210 cycle is populated (notice that the orbit goes very close to 
the edges of the simplex at later times). 



proaches the 310 face but eventually species 2, that was 
only apparently extinct, increases its density once again 
and the system stabilizes on the full 3210 heteroclinic cy- 
cle. 

To summarize, the system tends to decrease the amount 
of hierarchy (because of the x weighted connections) by 
the exponentially fast extinction of one species, converg- 
ing to a fully intransitive, non hierarchical, three species 
system [35j . Indeed, in mean field, any amount of tran- 
sitivity (measured by x) destroys the possible coexistence 
state that exists when x = 0. 

3. Simulations 

The dynamics on a lattice may be very different from the 
evolution predicted by the mean field equations, mainly 
because the range of interaction being much smaller than 
the system size, local correlations play an important role. 
Moreover, unless the system is very large, finite size effects 



exist and introduce stochastic effects. As an example, the 
invariants discussed in the previous section, quantities that 
are kept constant during the motion along closed orbits, no 
longer persist for finite systems, and density fluctuations 
eventually drive the system, through extinctions, into an 
absorbing state. These finite size effects become less im- 
portant for large systems and disappear for L — > oo, where 
L is the system linear size. In order to study the system 
on the lattice, we consider a square grid with N = L 2 
sites with periodic boundary conditions and, with the same 
probability, one of the four species is randomly assigned to 
each of those sites at t = 0. One site and one of its neigh- 
bors is chosen at random and the stronger site invades the 
other, depending on the species, with probability either 
unity or x- This step is repeated N times, what defines 
the time unit. Analogous to the mean field approach, the 
densities oscillate in time, however, the amplitude of these 
oscillations seems to decrease with the size of the system 
and tend to disappear for very large systems. 

Even though a deterministic system may have stable 
coexistence states, in its stochastic counterpart the fi- 
nite number of interacting agents induce fluctuations that, 
given enough time, eventually lead to the extinction of one 
or more species. However, as the system size increases, dis- 
tinct dynamical behaviors may be observed depending on 
the value of x- The dependence of the average character- 
istic time for an extinction to occur on the system size N 
allows for a classification of the possible occurrying scenar- 
ios [1, 55? ,H(|. The coexistence is said to be stable when 
the related deterministic dynamics presents a stable atrac- 
tor in the coexistence phase, and this is associated with 
an exponentially increasing time for the first extinction 
to occur as N increases. Analogously, the unstable state 
presents a logarithmic increase of the extinction time and 
the deterministic system approaches an absorbing state. In 
between, a power law dependence of the extinction time 
on the system size is related with the presence of closed, 
neutrally stable orbits in the deterministic case. The top 
panel of Fig. [4] shows, for small values of x an d several 
linear sizes L, the probability that the system does not 
suffer any extinction up to the time t, Po(t), that is, the 
probability of a persistent coexisting state. The larger the 
system is, the longer it takes for Pq to start dropping. We 
may define a characteristic time for the first extinction, 
t(N), as the time when Pg drops to half its initial value, 
that is, Po(t) = 1/2. In the inset of Fig. gj top panel, 
one can observe, for the range of sizes considered here, 
that t(N) has an exponential growth and even for modest 
sizes, the time of the first extinction is very large. Extinc- 
tion 54] in this case is driven by very rare fluctuations and 
the coexistence is said to be stable HIH? ,Hi|. On the 



other hand, for large x, inset of Fig. [H bottom panel, the 
extinction time growth is logarithmic in N and even for 
very large systems (one order of magnitude larger than in 
the previous case), r is rather small. Coexistence in this 
case is unstable and even small fluctuations are able to 
drive some species to extinction 5^, 56j |. Thus, comparing 
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Figure 4: (Top panel) Probability of all species surviving up to time 
t, for x = 0-31 < Xc and several linear sizes L. Data are averages over 
at least 200 samples. The characteristic extinction time t, defined 
as Po( r ) = 1/2, is shown in the inset and grows exponentially with 
the system size. For comparison, we also show in the inset the char- 
acteristic extinction time for a finite fully mixed system (see text). 
(Bottom panel) The same, but for \ = 0.5 > Xc The characteristic 
extinction time shown in the inset has a logarithmic dependence on 
N. Notice the very different horizontal scales in the two panels. 



these two cases, there must be a dynamical critical value 
of x, Xc, separating those two quite distinct dynamical be- 
haviors of t(N), and a rough estimate places this critical 
value at Xc — 0.35. Indeed, as shown in Figs. E] and O 
t(x < Xc) and t(x > Xc) have very distinct asymptotic 
behavior [57] . While for x < Xc the mean extinction time 
t grows exponentially, above Xc this growth is logarithmic 
in N. The intermediate region, for x — Xc, the scaling of 
t with the system size is polynomial. For x > Xc species 
2 goes extinct and the three remaining ones converge to 
densities close to the fixed point Eq. ([7]). It is also im- 
portant to stress that a second extinction, when it occurs, 
takes a much longer timescale. Both insets of Fig. @] also 
show the comparison with the correspondent r for a finite 
fully mixed system. To simulate such a system, in each 
MC step, new neighbors are randomly assigned to each 
site, without any distance constraint. For both x > Xc 
and X < Xc, r grows logarithmicly with L, but with a 



Figure 5: Average extinction time r versus N for several values of 
X- The critical value of x, Xc, is slightly above x = 0.35. Notice 
the two different asymptotic behaviors, r(x < Xc)/N — ¥ oo and 
T (x > Xc)/N — > 0. For x — Xc, T is linear in N. The lines are only 
guide to the eyes. 



small declivity, and no distinction exists between the two 
regions. 

For x < Xc, due to the exponential growth of r, the 
coexistence state is said to be stable and large systems stay 
in a state in which all four species attain a non zero fixed 
point. The average asymptotic density 'p i for large systems 
can be obtained by extrapolating the above behavior, that 
is, Pi = linij-^oo limjv-^oo Pi(t) (notice that the limits are 
not interchangeable). The results are shown in Fig. [H] as a 
function of x- The dynamical transition is clearly seen as 
the point at which species 2 goes extinct. Notice that all 
four densities arc different, both above and below Xc- 
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Figure 6: Stationary densities as a function of x- Although on a 
lattice we do not observe an oscillating behavior, there are sample 
to sample fluctuations and we average over at least 1000 samples 
(symbols). The lines are the fixed point, Eq. 0, of the mean field 
equations. Notice that albeit the reasonable agreement for X > Xc, 
below this value the system is in a state, not captured by the mean 
field approach, in which all four species coexist. 
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Comparing with the mean field predictions, we notice 
several fundamental differences. First of all, the coexis- 
tence state only appears on the regular lattice since for 
X 7^ 0, in mean field, species 2 always gets extinct. Thus, 
the dynamical transition that we observe does not exist 
for a fully mixed system. Indeed, we can simulate the sys- 
tem with annealed neighbors, and the characteristic time 
of the first extinction is small, presenting a logarithmic 
growth with iV for all x, as can be seen in both insets 
of Fig. |U The second difference is that although oscil- 
lations are observed for finite systems, they tend to dis- 
appear for very large sizes [33| . A possible explanation 
is that for a sufficiently large system, several different re- 
gions will evolve almost independently with uncorrelated 
phases, such that the overall system no longer presents 
oscillations. The third difference is that while the mean 
field behavior is monotonous on \ (p 2 is always zero, p 
is always increasing and ~p x = p 3 is always decreasing) , on 
the regular lattice the behavior is non monotonous, Fig. [5] 
Species 1 and 3 densities, now resolved, increase for x U P 
to % c and decrease afterwards (species has the opposite 
behavior). Notice, however, that although above \ c the 
agreement with mean field is reasonable, below Xc there 
is both qualitative and quantitative disagreement. In par- 
ticular, all four densities are different and no fixed point 
predicted by the mean field approach has such a property 
for X ^0. " 

4. Conclusions 

We studied a minimal model for a trophic network pre- 
senting multiple loops of interacting species, focusing on 
the effects of a tunable transitivity on the persistence of 
the coexistence state. As the invasion rate x is changed, 
we observed two distinct dynamical phases separated by a 
transition at Xc — 0.35, one for x < Xc in which the coex- 
istence state is stable (the mean extinction time exponen- 
tially grows with the size of the system) and the other for 
X > Xc in which one species goes extinct on logarithmic 
timescales and the system ends up performing a heteroge- 
neous RSP game. At the transition region between those 
two regimes, x — Xc, T presents a polynomial scaling with 
the size of the system. This transition, and the coexistence 
state observed in the simulations are not captured by the 
mean field approach. 

For x = 0, each species around the external four species 
loop has one prey and one predator. In addition, for x^O 
there are four internal loops with three species, two in- 
transitive (013 and 023) and two transitive (012 and 123). 
In this case, because of the even number of species, the 
number of predators and preys of each species may differ. 
Thus, depending on the arrows orientation, there are three 
possible choices for the number of preys (or, equivalently, 
predators): (3, 2, 1, 0), (3, 1, 1, 1) and (2, 2, 1, 1). We only 
considered the last structure, Fig. [TJ that is somewhat in- 
termediate between an intransitive and a hierarchical sys- 
tem. The larger is x, the less intransitive the system is 



and one can expect that the amount of coexistence will de- 
crease. However, we have shown that under the presence 
of spatial correlations and a not too large transitivity, this 
system may persist in a state of full diversity over expo- 
nentially large timescales. For larger levels of transitivity, 
on the other hand, the system eventually evolves into a 
three species hierarchical system, irrespective of the spa- 
tial structure [35j . We thus observe a dynamical transition 
between these two regimes on the spatially structured sys- 
tem, at Xc not captured by a mean field analysis. Notice 
that any extinction drives the system into an absorbing 
state and diversity, due to the absence of mutations, is an 
always decreasing quantity for this class of model. Below 
Xc, mean field is not a good approximation for the lattice 
dynamics of our system either quantitatively or qualita- 
tively. The threshold value of x also indicates that above 
a certain level of intransitivity (x > Xc), spatial correla- 
tions are no longer important and the system on a lattice 
is attracted to the mean field fixed points (the densities 
are non oscillating). Even though the pair approximation 
is assumed to be a better approximation than mean field, 
in our case it does not provide a better description of the 
dynamics, in terms of fixed points. For short times the 
dynamics of the PA does look similar to the mean field 
dynamics, since for large values of x the density of species 
2 drops to extremely low values, but for longer times the 
dynamics is dominated by a heteroclinic cycle involving all 
four species. 

Non monotonous responses driven by the spatial correla- 
tions are observed, while the mean field approach predicts 
a monotonous behavior as x changes. For x < Xc, species 
2 and 3 (and, analogously, and 1) respond in opposite 
ways: while ~p 3 increases with x, 7>2 decreases. The op- 
posite behavior was predicted in the mean field approach. 
On the other hand, above Xc the trends agree with the 
mean field prediction. Interestingly, in spite of present- 
ing opposite behavior when x increases, species and 1 
both become more aggressive. Since both predate on 2, 
this species has the smallest density (and becomes extinct 
in mean field). For x > Xcj P2 = ano - the remaining 
three species form a heterogeneous RSP game that obeys, 
both on the lattice and in MF, the usual "survival of the 
weakest" principle: as the invasion rate x °f the weakest 
species (1) increases, its density ~p 1 decreases, while the 
density of "the prey of the prey of the weakest" [l^|, in 
this case species 0, increases. For x < Xc since all four 
species survive, the "the prey of the prey of the weakest" 
principle (l9| must be modified because some species have 
multiple preys. Although species and 1 have a wider 
range of possible targets, they are less efficient since their 
overall success rate is less than 1 (1 + x < 2), and may 
be considered the weakest species. Species 2 and 3, on the 
other hand, fully overtake their preys. Nonetheless, the 
prey of the two weakest (species 2), itself stronger than 
them, goes extinct. Thus, although there is no obvious 
generalization of the above principle, one notice that the 
ambiguity in defining strong and weak in this case may be 
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raised by allowing all six parameters to be different, what 
may, in turn, allow for such an statement. It is also clear 
that statements like this will become more intricate as the 
number of species increases. 

Further questions arise for such systems. For example, 
in order to better understand the effects of different levels 
of transitivity, in particular to probe anomalous responses 
as the "survival of the weakest" , the study of other trophic 
structures with four species, and larger values of S as well, 
is important. In addition, finite populations may have a 
different behavior. Indeed, the community size, besides 
setting the scale for the average extinction time, may also 



influence which is the surviving species 58|. One may 
also probe the robustness of the results presented here, 
for example, by studying different lattices (random graph, 
small world, etc), dimensions and initial conditions. On 
a regular lattice, geometric and dynamical properties of 
the evolving groups are also of interest [3(1, [37| ■ How to 
properly quantify the transitivity of a trophic network and 
correlate it with the coexistence present in a population is 



still an open problem (e.g., [32l. l33l. l59H62l | and references 
therein). We have shown that the structure alone is not 
enough to predict whether there will be coexistence or not. 
Considering the trophic relations as a weighted network 
may lead to an index allowing different levels of coexistence 
based on the same structure. Finally, allowing general 
weights on the trophic network [r| shall present an even 
richer behavior in the presence of crossed interactions. 
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